Optimal superluminal systems 
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We demonstrate that significant effects in the "superluminal propagation" of light-pulses cannot be 
observed without involving systems whose gain explodes outside the pulse spectrum. We explicitly 
determine the minimum norm of the gain to attain given superluminal effects and the transfer 
function of the corresponding optimal system. The gain-norms which would be required with the 
most efficient systems considered up to now (dispersive media, photonic barriers) to attain the same 
effects are shown to exceed the minimum by several orders of magnitude. We finally estimate the 
largest superluminal advances which could be attained in a realistic experiment. 
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PACS numbers: 42.25.Bs, 02.60.Gf, 03.65.Xp 

The apparently superluminal propagation of light- 
pulses has been observed with various systems, mainly 
with systems involving anomalously dispersive media or 
photonic barriers. For reviews, see, e.g., Q, Q> IH E 01 • 
In these experiments, the envelope of the pulse having 
covered some distance L is nearly identical to that of the 
incident pulse and in advance of that of a pulse which has 
covered the same distance L at the velocity c of light in 
vacuum. This surprising behaviour is not at odds with 
the relativistic causality. Indeed the signal received at 
some time t is not the consequence of the signal emitted 
at a well-defined time but of all the signals anterior to 
t by more than L/c. Otherwise said, there is no cause- 
to-effect relation between the homologous points of the 
envelopes of the incident and transmitted pulses and the 
widespread statement that the pulse maximum leaves the 
system before it even enters it is somewhat misleading. 
The phenomenon is however quite puzzling and keeps the 
subject of an intense theoretical and experimental activ- 
ity. 

In fact Mother Nature resists to a violation of her 
principles even when this violation is only apparent and 
convincing experiments of superluminal transmission are 
very difficult to achieve. By convincing experiments, we 
mean experiments where (i) the envelopes or the intensity 
profiles of the pulses arc detected in real-time and true- 
shape (ii) the transmitted pulse is directly compared to 
the pulse having propagated at the velocity c (iii) the 
superluminal advance a is large compared to the opti- 
cal period (iv) it is significant with respect to the pulse 
duration, say larger than 10% of the full width at half 
maximum (FWHM) of the intensity profile r p (v) the 
pulse distortion (including noise and parasitic signals) 
is small compared to the relative advance a/r p . Note 
that (iii) is a consequence of (i) since the real-time de- 
tection of the envelope requires a time-constant large 
compared to the optical period. There are few exper- 
iments meeting, even approximately, the previous con- 
ditions HH&HHIlllElEl- Though all-optical 



experiments are possible, only hybrid systems have been 
used up to now. They combine an optical part, respon- 
sible for the superluminal effects, and a wide-band elec- 
tronic device whose function is to normalise the ampli- 
tude of the transmitted pulse. In most experiments, the 
transmission of the optical part, usually a resonantly ab- 
sorbing medium 0, 0, 0, El 0] or a photonic barrier 
0- [HI , is low and the electronic device is an ampli- 
fier. To our knowledge, only one experiment [13] has 
evidenced significant superluminal effects with an active 
optical part (amplifying medium). The normalisation is 
then achieved by a suitable attenuation. In the following, 
we naturally include the normalisation device (amplifier 
or attenuator) in the system under consideration. 

As already noted in previous papers dealing with par- 
ticular arrangements (see, e.g., p|), large superluminal 
effects are only attained with systems whose gain ex- 
plodes outside the pulse spectrum. We will show that this 
is true for any physically realisable system and determine 
the lower limit to the gain norm required to observe given 
superluminal effects. This result is of special importance 
since in a real experiment the gain-norm should be lim- 
ited to avoid problems of noise (no matter its origin), of 
instability and of hypersensitivity to parasitic signal and 
to localised defects in the incident pulse profile [2j . Con- 
versely restricting the gain to realistic values determines 
the upper limit to the actually observable effects. 

The problem is studied in the frame of the linear sys- 
tems theory [l^. We denote by e(t) and s(t) the en- 
velopes of the incident and transmitted pulses and by 
E(u>) — J_ e(t) exp(—iujt)dt and S(u>) their Fourier 
transforms. The envelopes are assumed to be slowly 
varying at the scale of the optical period. Their Fourier 
transforms are then concentrated around in a region 
of width small compared to the optical frequency. In 
all the sequel, t designates the local time, equal to the 
real time in e(t) and retarded by the luminal transit time 
L/c in s(t). The system is characterised by its impulse 
response h(t) or its transfer function H(u>), such that 
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s(t) = h(t) ® e(t) and = H{uj)E{lo). We assume 

that and H(uj) have a finite energy and that H(u>), 

Fourier transform of h(t), has a continuation H(z) in 
the complex plane [z — 



— 



pe ). In our lo- 
cal time picture, the relativistic causality imposes that 
h(t < 0) = 0. Otherwise said, H{z) belongs to L 2 (R), 
the Hilbert space of functions F(z) square summable 
on the real line R endowed with the norm \\F\\ R such 
that ||F|| R = \F(uj)\ 2 duj and, more precisely, to the 
Hardy space H 2 ) of functions F analytic in the lower 
half-plane II_ (y < 0) which are Fourier transform of 
some causal function / £ £ 2 (0, oo) [lfij . 

We want s(t) to be as close as possible to e(t+a) where 
a is the superluminal advance (a > 0). In L 2 norm, the 
distortion is defined by 



D 



where H„ 



\e(t- 



|e(i)|| 
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(1) 



lula is the transfer function of the non causal 
system perfectly achieving the advance a without any 
distortion. With a real (causal) system, the distortion 
will be low if H(uS) w H a (uj) in the region around u> = 
where |_E7(o7) | is concentrated. 

To keep tractable calculations, we consider the case 
E(uj) — Eq for \uj\ < uj c and elsewhere. By taking E = 
7r and lj c = 1, this amounts to take as reference a pulse 
of intensity profile \e{t)\ 2 = sin 2 t/t 2 (FWHM t p = 2.78). 
The distortion then reads D = \\H a — H\\j /y/2 where 
\\F\\j denotes the norm L 2 of F restricted to I = [—1, 1]. 
In the situations of physical interest D<1 and ||i?|| R = 
\\Hf I +\\H\\ 2 T ^2+\\H\\ 2 J wheveJ = [-oo, -1]U[1, oo]. In 
this model, the problem may then be stated : given a > 
and D > 0, minimise Q 
H £ H 2 (n_) and \\H a - H\\ T < L>V2 

Based upon a conformal map that sends the unit disk 
B (p — 1) onto the lower half-plane, we introduce the 
map F = ^ (F) defined by 



if II j under the constraints 



F(z) = viz (F) (z) = 



F 



l-z \2i{l-z) 



(2) 



It is an isometry from L 2 (M.) to the Hilbert space L 2 (T) 
of the unit circle T endowed with the norm ||.F|| T such 

that \\F\\j — J 27V \F(e l9 )\ 2 d9/2ir. It sends the subspace 
_ff 2 (n_) onto the corresponding Hardy space 7J 2 (B) of 
the unit disk B . We denote by I and J the subarcs 
of T, transforms of I and J by the map ^ . Then this 
map allows one to restate the problem in the unit disk B 
instead of the lower half-plane : given a > and D > 0, 



minimise Q = 

H a -H 



H 



ano 



under the constraints H £ H 2 (B) 



< Dt/2. 
i _ 

Stated with a general function K £ L 2 (I) instead of 
the particular H ai this question has been originally con- 
sidered in [lr| and more recently in with important 



extensions. The solution H opt of the problem exists and is 
unique. Note that, in our case (K = H a ), the constraint 
\\H a - H\\j < Dy/2 is saturated, i.e. \\H a - H\\ z = Dy/2. 
The solution H op t can formally be written under the an- 
alytic form [17|: 
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Opt 



(l + \$)- l P H2 (H a ) 



(3) 



In this expression H a is defined as H a on I and on J, 
Pfj2 denotes the orthogonal projection from L 2 (T) onto 
H 2 {p) and $ is the so-called Toeplitz operator 17] acting 
on H 2 (B). It is such that = P H ^(F) where F is 

defined as F on J and on I. Finally A £ [— l,oo] is 
an implicit parameter. It is the unique real number such 
that \\H a - H\\j = Dy/2. 

From a computational viewpoint, i t ap pears natural to 
consider Q and D as functions of A [17J. It follows from 
Eq|3]that Q and D respectively increases and decreases 
as A decreases. As A — > —1, Q — > oo while D — > 0. In 
physical terms, this confirms that a low distortion will 
always be paid at the price of a large gain-norm. We 



have then \\H, 
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Given a and D, the previous analysis leads to the fol- 
lowing algorithm for the computation of the minimum 
gain norm Q and the corresponding function H opt : (i) 
Choose A < — 1 and compute H opt given by EqEl (ii) 
Compute D. If it is too large (resp. small), decrease 
(resp. increase) A. Goto(i). Such a dichotomy algorithm 
has been implemented in the software package Hyperion 
developed at INRIA (Institut National de Recherche en 
Informatique et Automatique) by the APICS team [Isj ]. 
See also |lj| for a closely related algorithm. EqEl which 
is infinite dimensional, is approached by truncating the 
expansions of the involved functions so as to consider 
only their Fourier coefficients of indices —N < j < N. 
The optimal transfer function H opt (tu) is finally obtained 
by inverting Eq|2 



H opt (uj) 



2iu + 1 " 



-H, 



Opt 



1 2iuj + 1 



(4) 



Note that H opt (uj) behaves as 1/iui for \uj\ — > oo. This 
behaviour is that of a first order filter as used in every 
detection chain. Any further filtering of the high fre- 
quencies will obviously damage the performances of the 
system. To close this short presentation of our minimi- 
sation procedure, we remark that it mainly lies on the 
separation of the spectral domains where the distortion 
and the gain-norm are computed. We have chosen the 
pulse profile leading to the simplest calculations but the 
procedure might be adapted to any pulse provided that 
its Fourier transform has a compact support. 

Calculations of the minimum gain- norm Q, of the cor- 
responding transfer function H opt (uj) and of the transmit- 
ted signal s(t) have been made for a/r p (resp. D) ranging 
from 0.36 to 2.2 (resp. 2 to 30%). Satisfactorily enough, 
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FIG. 1: Amplitude-gain G and phase ip (radian) of the opti- 
mal system as functions of the frequency. Parameters: a = t p 
and D = 15%. 
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FIG. 2: Frequency-dependence of AG = G — 1 and of A<p = 
p — una in the useful band. The group advance a g deviates 
from a by Ao = d{Ap)/duo \ u =o, that is Act ~ —0.40 and 
Aa/a ~ —14%. Parameters as in Fig0 



the optimal system would allow one to conciliate signif- 
icant advance, moderate distortion and reasonable gain. 
For instance a = r p with D = 15% would be obtained for 
Q = 100. Fig^ shows the overall frequency-dependence 
of the amplitude-gain G (uS) = \H opt {io)\ and of the phase 
ip (lu) = arg [Hopt(u))] in this reference case. As expected, 
the gain reaches its peak-value near the frontiers of the 
"stop band" (in fact the useful band for superluminal 
systems). The short ringing close to these frontiers orig- 
inates from the finite number of Fourier coefficients used 
in the calculations (N = 2000). The asymptotic values of 
the phase are <p — ±9n/2 for lu = ^foo, in agreement with 
EqQ] The extra phase-rotation of 87r entails that H opt (z) 
has four zeros in the half-plane y < and, consequently, 
that H opt {uj) is not minimum-phase [ljj. The differences 
AG = G - 1 and Atp = tp - uua for -1 < lu < 1 (Fig© 
illustrate how H opt {oj) deviates from the ideal transfer 
function H a = e luja in the useful band. We remark that 
the group advance a g = dtp/cLu | w= o differs from the ef- 
fective advance a by an amount approximately equal to 
the distortion (in our local time picture a g = L/c— L jv g 
where v g is the group velocity). Finally, the envelope 
s(t), inverse Fourier transform of H opt (ui)E(uj), and the 
intensity profile |s(i)| 2 of the transmitted pulse are dis- 
played FigEI 

The efficiency of a superluminal system may be char- 
acterised by its ability to achieve given effects with gains 
as small as possible. As above-noticed, the gain of all 
the optimal systems has the same asymptotic behaviour 
(G cc 1/lo) and reaches its peak- value M near u> = ±1. 
Consequently Q and M are roughly proportional and can 
indifferently characterise the system gain. The peak-gain 
M, independent of the frequency scaling, is retained in 
the sequel. This choice facilitates the comparison of the 
optimal systems with the most efficient systems used or 
proposed up to now. Since high optical gains exaggerate 
the problems of instability and noise (amplified sponta- 
neous emission) and are difficult to achieve with the suit- 
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FIG. 3: Intensity profile of the pulse transmitted by the op- 
timal system (full line). The profiles of the incident pulse 
advanced by a — r v (dotted line) and the main lobe of the 
incident pulse (dashed line) are given for reference. Insert: 
Envelopes s(t) (full line) and e(t + a) (dotted line). Parame- 
ters as in Fig0 



able frequency-profile [13| , we restrict ourselves to sys- 
tems whose optical element, responsible for the superlu- 
minal effects, is passive. More specifically, we consider a 
dilute medium involving (a) an isolated absorption-line 
H H ES El E3 or (b) a doublet of absorption-line [| 
and an uniform Bragg-grating written (c) on a classi- 
cal optical fibre or (d) on a hollow fibre. Since all 
these elements are almost transparent outside the low- 
transmission region (the useful band), the peak-gain M 
is nothing but the gain of the amplifier used to normalise 
the amplitude of the transmitted pulse. The transfer 
functions are optimised by adjusting the system param- 
eters with a genetic algorithm. A rapid convergence is 
obtained by starting the calculations with initial values 
such that H(0) = 1 and a g — a. For the doublet (b), 
the initial value of the line-splitting is chosen such that 
the 2 nd order distortion cancels 0. Fig0| shows the re- 
sults obtained for a reference distortion D = 15% and 
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M ranging from 10 to 3 x 10 4 . No need of a lens to see 
that the optimal system is much more efficient that the 
systems (a), (b), (c) and (d) to attain large superluminal 
advances. For instance, a peak-gain M = 84 theoretically 
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FIG. 4: Relation between the peak-gain M and the relative 
advance a/r p for a given distortion (D = 15%). (o) stands 
for the optimal system while (a), (b), (c) and (d) respectively 
relate to the so designated systems (see text). 

suffices to observe an advance a = t p with D = 15% (Fig 
HJ but values as large as 1600, 3400, 6.4 x 10 6 and 4.9 x 10 7 
would be required with the systems (b) , (d) , (a) and (c) 
respectively |2j}. The latter dramatically increase if a 
lower distortion is required. Again for a = t p but with 
D = 7% they raise to 7.9 x 10 4 , 2.1 x 10 7 , 2.3 x 10 14 
and 4.9 x 10 15 while M only reaches 174 for the optimal 
system. By comparison, we stress that achieving experi- 
ments with systems whose peak amplitude-gain exceeds 
10 4 is absolutely unrealistic. 

The situation is much less catastrophic when one ex- 
amines the superluminal effects which can be attained 
for a fixed peak-gain. Taking M — 1000 (realisable in 
a careful experiment) and D = 15% as reference val- 
ues, Fig0] shows that the relative advance a/r p attained 
with the simplest arrangement (medium with an isolated 
absorption-line) is only 2.4 times below the theoretical 
limit (1.6) and that the ratio falls to 1.7 by involving 
a line-doublet. Using non uniform fibre-Bragg-gratings 
could further reduce this ratio. Indeed, at least in prin- 
ciple, these elements allow one to synthesise any transfer 
function in transmission as long as it is minimum-phase 
[23. This restriction entails that the optimal transfer 
function (not minimum-phase) and thus the upper limit 
to the advance could be approached but not equalled 
with these systems. The same remark applies to the dis- 
persive media whose transfer function is the exponential 
of a causal function and is thus also minimum-phase 
Anyway, whatever the system is, superluminal advances 
exceeding two times the full width at half maximum of 
the pulse intensity-profile are unattainable. 
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